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^ . ABSTRACT 

00 ' Context. About two dozen exo-solar debris systems have been spatially resolved. These debris discs commonly display a variety of structural 
T-H , features such as clumps, rings, belts, excentric distributions and spiral patterns. In most cases, these features are believed to be formed, shaped and 

. maintained by the dynamical influence of planets orbiting the host stars. In very few cases has the presence of the dynamically important planet(s) 

1 I • been inferred from direct observation. 

P^ ■ Aims. The solar-type star q' Eri is known to be surrounded by debris, extended on scales of i 30". The star is also known to host at least one planet, 
^O , albeit on an orbit far too small to make it responsible for structures at distances of tens to hundreds of AU. The aim of the present investigation 
is twofold: to determine the optical and material properties of the debris and to infer the spatial distribution of the dust, which may hint at the 
presence of additional planets. 

Metiiods. The photodetector array camera and spectrometer (PACS) aboard the HerscheJ Space Observatory allows imaging observations in the 
/i>) , far infrared at unprecedented resolution, i.e. at better than 6" to 12" over the wavelength range of 60yum to 210yum. Together with the results from 
5^ 1 ground-based observations, these spatially resolved data can be modelled to determine the nature of the debris and its evolution more reliably than 
"*T* ■ what would be possible from unresolved data alone. 

jrt " Results. For the first time has the q' Eri disc been resolved at far infrared wavelengths. The PACS observations at 70 /jm, 100 yum and 160 /jm 
' ' , reveal an oval image showing a disc-like structure in all bands, the size of which increases with wavelength. Assuming a circular shape yields the 
. inclination of its equatorial plane with respect to that of the sky, / > 53°. The results of image de-convolution indicate that i likely is larger than 
T-H • 63°, where 90° corresponds to an edge-on disc. 
^ ' Conclusions. The observed emission is thermal and optically thin. The resolved data are consistent with debris at temperatures below 30 K at radii 
C , larger than 120 AU. From image de-convolution, we find that q' Eri is surrounded by an about 40 AU wide ring at the radial distance of ~ 85 AU. 
CO . This is the first real Edgeworth-Kuiper Belt analogue ever observed. 
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(<— ^ ' 1. Introduction implying that the whole sample covers a decade in stellar mass 

— I ' , from 0.2 to 2 Mq. 

We report new observations of the debris disc around q Eri _,, „.„„ it: — :: — ; , i ^„.^ , . ,„„ 

/iir-.in/c/n iir.cn/c\ .u it ; ; o <-.u . The PACS flPogutsch et alJ2010l) observations at lOOum aim 

(HD 10647, HR506) using the Herschel Space Observatory , , . '— r^ r; f ", , , j- ■ 

l r>iu J ..1 l l oninK t-u u .• c . f i r^ at the detection of the stellar photospheres down to the contusion 

vj JPilbratt et al.ll201(JI) . The observations form part of a larger Kev . ., . , • • ;o,xrx ,■ , ^ rr, 

^r' ■ y, ,T^-n\ ■ T>TTMT-cfn uujuj- J noisc With a signal to noise ratio (S/N) of at least 5. Together 

rN . Programme (KP), viz. DUNE3T which IS described in more de- ., , . ° . , , ,, , ,/ , ,. .,, , ■ 

$-H V .,? I t^- ^ 1 N -,rvi/4 . TT • 1. • r . . withobservationsintheother/Zeracne/bands, this will lead to an 

rrt I tail bv lEiroa et al.l (120 101) . Here, we give a brief summarv to put , , , ^ ,. ,.,,,, ■ ., ■ 

'•^ ,, ~ — I — TTu^~~i — n • . T t T-u i^iT^Ti-c r^n • Unprecedented charactcrisation ot discs and Will allow detailed 

the contents of this Letter into context. The DUNES KP is a , ^ . , , „. , . ^ 

.,. -^ I- •. J » J -^u ^u ^ c A- ■ J u theoretical modelling. It is foreseen that fractional luminosities 

sensitivity hmited study with the goal of discovering and char- . ,. ^ ^ .° ,„_7 .,, , , , . . ., , 

^ . . ^ , , ? »u T-j .u T> • n 1. Ldust Lnj of a few times 10 will be reached, i.e. similar to that 

acterismg extra-solar analogues of the Edgeworth-Kuiper Belt ""?_' ^ r ,u c ^ c . ,q. inn<; i » > i I hnnnK 

fTDVD\ ■ u- A ^ ^- »■ 1 1 ? I. TD ^ A u of thc EKB of thc Solar Svstem Stern 1996; Jewitt et alJI2009l) 

(EKB) in an unbiased, statistical sample of nearby F, G and K , , , ^ ■ — \ — '■ Itt — 

rru 11 1- V J -.u A- and a more than order of magnitude improvement over Spitzer 

main-sequence stars. The sample is volume limited, with dis- , 17; — -. nrrrrr^ & r i- 

< or, A u A f , u f n 1 data dBrvden et al.ll2009h . 

tances S 20 pc, and spans a broad range of stellar ages, from ~ 0. 1 — 

to roughly 10 Gyr. In addition to the object of the present study ^ main-sequence star with an infrared excess larger by more 

(q' Eri), a number of M- and A-type sta rs will be observed i n *an three orders of magnitude than this limit is the late F- 

coUaboration with the DEBRIS-KP team (iMatthews et al.l201flh . type star q En. This will potentially allow to study the mate- 

rial giving rise to the excess in great detail. The star is known 

* Herschel is an ESA space observatory with science instruments to be accompanied by^_giant planet, q' Erib, orbiting at the 

provided by European-led Principal Investigator consortia and with im- distance of 2 AU (.Butler et al. 2006, and references therein), 

portant participation from NASA. which corresponds to OV 1 at the distance of the system (17.35 pc; 



DUst around NEarby Stars, [http://www.mpia-hd.mpg.de/DUNES/| iPerrvman et alJfT997t Ivan Leeuwenll2007h . Such a small angle 
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Fig. 1. PACS photometric imaging of q' Eri at, from left to right, 70 //m (blue), 100 yum (green) and 160 yum (red). The 70 //m image 
was taken in chop-nod mode, whereas the other two in scan-map mode. The upper panels display the reduced observations. Below, 
de-convolved images are shown, using o bservations of a Boo for the definition of the PSF. Displayed are the results for ten iterations 
of a MEM algorithm (iHolUs et aljl992l) . The star defines the origin of the frames, i.e. offset coordinates (0, 0). Within the positional 
accuracy (2" rms), the stellar position and the centre of the elliptical brightness distributions coincide (see Table[T]) and offsets are in 
seconds of arc. The lowest contours are at 5% of the maximum values and consecutive steps are also by this amount. At the distance 
of the star, 20" corresponds to 350 AU. 



will not be resolved by the observations described in this pa- 
per, but q' Eri is also surrounded by a ring or belt system of 
scattering and thermally emitting dust particles. The extent of 
the system, which from optical to sub-millim etre wavelengths 
is up to several tens of ar c seconds in size (e.g.. lStapelfeldt et al.l 
l2007tlLiseau et al.l2008l and references therein), should be read- 
ily accessible to PACS at the Herschel telescope (flj^jj - 7 x 
/l^m/100//m). 

This contribution presents the observed properties of the 
q' Eri system by Herschel. The results of theoretical model cal- 
culations will be communicated by Augereau et al. (in prep.). 
The observations and data reduction are described in Sect. 2, 
with the results presented in Sect. 3. These are discussed in 
Sect. 4 and, finally in Sect. 5, our main conclusions are briefly 
summarised. 



2. Observations and data reduction 

During this initial observing run, two different modes of ob- 
serving were executed for test reasons, in order to optimise the 
efliciency of the programme in terms of observing time and 



signal-to-noise ratio (S/N). These modes were, respectively, the 
chop-and-nod mode (70, 100 and 160 jum), adopted for point 
source observing, and the scan-map option (100 and 160//m), 
for extended sou rces (Fig.[T]i. These modes are described by 
lEiroa et al] (1201 Ol) and in a future technical note. 

The reduction was done within the Herschel interactive pro- 
cessing environment, using HIPE_v2.0.0_RC3, and with scripts 
for the pipeline developed by members of the PACS-ICCQ- At 
the medium scan speed of 20" s"', two maps along position 
angles 63° and 117°, respectively, were obtained in order to 
minimise striping in the resultant images. At both 100 //m and 
160jum, the sky noise was lower in the scan-map data (2.8 and 
5.2) X 10"^ Jy arcsec"^, respectively, as compared to (3.9 and 
12) X 10"^ Jy arcsec"" of the chop-nod observations. These val- 
ues are comparable to the surface brightness of the scattering 
disc at 0.6 //m (1-3 x 10~^ Jy arcsec"-, as deduced from data in 
IStanelfeldt et"al]l2007h . 



^ Herschel science demonstration phase data processing workshop, 
14-16 December 2009, ESAC, Spain. 
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Table 1. Brightness, orientation and extent of the q' Eri debris system 



Observation 
ID and Mode 



ARA, ADec" 



Flux'' 
(mJy) 



2a'' 
(") 



2b'' 



PA' 

D 



n 



Alternate flux 
measurement (mJy) 



1342187142 - 70 : chop-nod 

1342187141 - 100 : chop-nod 
134218739/40 - 100 : scan map 

1342187142 - 160 : chop-nod 
134218739/40 - 160 : scan map 



0.9, -0.9 


828 ± 83 


13.9(11.0) 


7.9 (4.4) 


54.4 (60.0) 


55.4 (66.4) 


-1-2,0 


816 + 82 


15.5 (12.6) 


8.8(5.1) 


54.4 (56.2) 


55.4(66.1) 


0, -1.4 


810 + 81 


15.1 (12.5) 


8.8 (4.7) 


56.1 (54.8) 


54.4 (67.9) 


1.4,-1-0.9 


529 + 106 


18.9(13.5) 


12.3 (7.8) 


51.3(48.9) 


49.4 (54.7) 


1.4,-1.6 


537 + 107 


19.3 (13.3) 


12.5 (6.5) 


51.2(55.7) 


49.6 (60.7) 



859+ 6 Spitzer/MlPS'- 
1080 + 36 IRAS FSC 
1080 + 36 IRAS FSC 
453 + 50 Spitzer/MIPSf 
453 + 50 Spitzer/MIPS^ 



Notes to the Table: 

° Offsets of ellipse centre relative to stellar coordinates: proper motion corrected to IRCS 2009.89 (+0.186s and - 1'.'055 iPerrvman et al.ll 19971 : 

Ivan Leeuwei3l2007h 0l''42™29»502, -53° 44' 28'.'06. The photospheric fluxes are I7mJy at 70/jm, BmJy at 100 yum and 3mJy at 160 yum (cf 

Fig.|2j. 

* Flux for Gaussian ellipse with fitted FWHM to major and minor axis, 2a and 2b, respectively. Values in parentheses refer to de-convolved images. 

'^ Position angle measured from North over East. 

'' Lower limit to the i nclination, where ; = 90° refers to an edge-on geometry. 

'' lTrillinget"anil2008h . 

^ iTanneret'aiT ( l2009l) . 



The prime calibrator Arcturus (a Boo; e.g., ICohen et alJ 

l2005h was observed close in time with q' Eri and was used to 
provide the instrumental point spread function (PSFJj. 

The filters at the reference wavelengths of 70, 100 and 
160 fim are referred to as blue, green and red, spanning 
60 - 85 ^m, 85 - 1 30yum and 130 - 210jum, respectively 
(iPoglitsch et al.ll2008h . The PSF has a tripolar shape at low in- 
tensities, i.e. a few percent of the peak value, and the half-power- 
width is circular for the 70 //m and 100//m filters, but somewhat 
elongated at 160/Ym in the scan direction. Currently, the esti- 
mated accuracy of the absolute flux calibration is estimated to 
be better than 20% in the long- wave and better than 10% in the 
short-wave bands. 



3. Results 

In the Scan Map AOT release is the following note: The fluxes 
are too high and have to be scaled down by the following factors 
1.05 in the blue band, 1.09 in the green band 1.29 in the red 
band, and in Table[Tl no colour correction has been applied to 
the reported fluxes. The internal consistency of the PACS data 
- at a given wavelength but for different observing modes - is 
strikingly good and lends confidence to the quality of these data. 
From the comparison with previous measurements, it is apparent 
that IRAS fluxes are on the high side, whereas long-wave Spitzer 
data are somewhat on the low side. 

At the level of the pointing accuracy of Herschel (2" rms, 
corresponding to 35 AU), there is no significant offset of the cen- 
tres of the elliptical isophotes with respect to the position of the 
star (Table[T]). This uncertaint y is much larger th an the offset of 
8 AU observed for Fomalhaut dKalas et al.ll2008h . For the nomi- 
nal wavelengths of 70jum, 100 jum and 160 jum, the angular reso- 
lution is limited to about 6", 7" and 1 1'.'5, respectively. The im- 
age of q' Eri seems clearly resolved, therefore, only in the long 
dimension. 

The average of the position angle of the elliptical source is 
PA - 54° and the average of the lower limit to the inclination 



^ See http : //Herschel . esac . esa . int/AOTsReleaseStatus . shtml 
http : //Herschel . esac . esa . int/Docs/AOTsReleaseStatus/PACS_Scan 
ChopNod_ReleaseNote.22Feb2818.pd£ 

http : //Herschel . esac . esa . int/Docs/AOTsReleaseStatus/PACS_Phot 
Map_ReleaseNote.23Feb2S18.pdf 
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Fig. 2. This figure is similar to the SED of q' Eri shown in 
iLiseau et all (l2008h with the difference that the PACS data of 
Table[T]are shown by the filled symbols. In addition, the 160 //m 
flux of Spitzer (Tanner et al. 2009j) has also been added. As be- 
fore, the solid curve for the IR excess is that of a single temper- 
ature blackbody (J3 = 0, where Ky oc v^), and the dashed line 
refers to the ring-belt composite model of that paper. 



of an intrinsically circular feature (ring or belt) is / > 53°, with 
estimated uncertainties of ~ ±5° (/ - 90° for an edge-on geome- 
try). For the de-convolved images (see below), the disc would be 
seen more edge-on and in this case, / > 63°. This is consistent 
with the tilt deriv ed from optical images of the light scattered off 
the disc (/ = 76°. IStapelfeldt et al.ll2007h . 
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-20 

Offset (arcsec) 

Fig. 3. One-dimensional cuts (averages of 5 pxl wide strips), 
along the major axis, through the 100 //m scan map image, from 
which the stellar source has been subtracted prior to the de- 
convolution. The histogram depicts the observed data, whereas 
the smooth line shows the result of the applied de-convolution 
algorithm. 



ysed sub-frame, the flux was conserved within 7% by the MEM 
routine. Prior to the de-convolution, a stellar point source with 
photospheric flux of 8mJy at IQQfim (Table[T]) was subtracted 
from the PACS image. The resulting sharpened image reveals a 
central broadened plateau and a central depression with a depth 
of about 2%, which is consistent with the debris residing in a 
ring or belt around the star. 

With standard assumptions regarding the emitting grains (as- 
tronomical silicates, a blow-out size limit flmin - 0.6 jum for the 
F9 V star and a -3.5 power law index for the size distribution) 
we find that the 100 jum surface brightness profiles along the ma- 
jor and minor axes are well reproduced assuming a disk incli- 
nation of about 70° and a two-parameter model for the surface 
density, £(r). These parameters are the peak density position, 
'"max, and the power law index of the surface density profile for 
r > '"max- The best fit to the surface brightness profiles is con- 
sistent with a ring-like disc, having values of r,nax ~ 85 AU and 
S(r > Tniax) 0^ '" "'i respectively. A more elaborate model, with 
the size distribution computed self-consistently and taking into 
account the profiles also at other wavelengths, will be presented 
by Augereau et al. (in prep.). 

This roughly 40 AU wide ring or belt at about 85 AU from 
the star appears similar to the EKB of the Solar System. Based 
on an analogy with the debris disc around Fomalhaut and on 
theoretical expectations, it is quite possible that another gas giant 
planet, q' Eric, could be orbiting the star inside the inner belt 
edge. Given the age of the system, i 2 Gyr, the direct detection of 
q' Eri c, for instance by means of coronogr aphy, can be expected 
to be hard (see, e.g.. lBeichman et al1l2006l) . 



4. Discussion 

4. 1 . The spectral energy distribution (SED) revisited 

Figure|2]is taken from the paper bv lLiseau et al.l (|2008|) . but with 
the PACS data included. The 70/Ym and lOOyL/m fluxes fall essen- 
tially on top of the 60 K blackbody curv e, shown as a solid line. 
For 160yum, the Spitzer observations of iTanner et al.l (l2009h are 
also included. These seemed to confirm the ring-belt model in 
our earlier paper The new PACS datum is marginally below the 
blackbody, but also marginally above the composite SED and, at 
the moment, we need to leave this issue as undecided. 

4.2. Image de-convolution 

From high-S/N data, possible fine-structural details in the im- 
age can be retrieved. This would require image sharpening tech- 
niques, such as, e.g., image de-convolution with the known PSF. 

The images were de-convolved with the rotated Arcturus- 
PSF using the algorithm for a max imum entropy method (MEM) 
developed by iHollis et al] ( Il992h . Guided by experience, we 
stopped the routine after about a dozen iterations in order not 
to produce artefacts and/or false super-resolution. More auto- 
matic and/or quantitative methods were tried, but were eventu- 
ally abandoned due to insufficient reliability. In Fig.[Tl the direct 
observations are compared with the de-convolved images and in 
Table[Tl values measured on the sharpened images are reported. 

In Figure[3] one-dimensional cuts along the major axis of 
the observed and de-convolved lOOyum image are displayed. It 
would have been more natural to select the observation at 70/im, 
having the highest spatial resolution, but the 100/im scan map 
data are of considerably higher S/N, outweighing the apparent 
resolution advantage of the shorter wavelength data. In the anal- 



5. Conclusions 

Based on imaging observations with PACS in the three photo- 
metric bands at IQfim, 100//m and 160yum we find that 

• the debris around the solar- type star q' Eri has an oval- 
shaped brightness distribution, the size of which increases 
with the wavelength; 

• the integrated flux density at these wavelengths leads to an 
SED which is in good agreement with earlier results; 

• the very high signal-to-noise of the 100 //m scan map is ad- 
equate to sharpen the image using an image de-convolution 
technique, revealing a ring-like structure with maximum sur- 
face density at ~ 85 AU from the star; 

• with a width of about 35 to 45 AU, this ring or belt around 
the F9 V star q' Eri is similar to the Edgeworth-Kuiper Belt 
around the Sun. This may hint at the presence of another 
planet, q' Eric. 
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